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Preamble 

As in earlier reports, we will continue to break our effort into six distinct units: 

• Atmospheric Correction Algorithm Development 

• Whitecap Correction Algorithm 

• In-water Radiance Distribution 

• Residual Instrument Polarization 

• Pre-launch/Post-launch Atmospheric Correction Validation 

• Detached Coccolith Algorithm and Post-launch Studies 

This separation has been logical thus far; however, as launch of AM-1 approaches, it must be rec- 
ognized that many of these activities will shift emphasis from algorithm development to validation. 
For example, the second, third, and fifth bullets will become alm ost totally validation-focussed 
activities in the post-launch era, providing the core of our experimental validation effort. Work 
under the first bullet will continue into the post-launch time frame, but will be driven in part by 
algorithm deficiencies revealed as a result of validation activities. We will continue to use this 
format for CY98. 
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Abstract 


Significant accomplishments made during the present reporting period: 

• We expanded our “spectral-matching” algorithm (SMA), for identifying 
the presence of absorbing aerosols and simultaneously performing atmo- 
spheric correction and derivation of the ocean’s bio-optical parameters, to 
the point where it could be added as a subroutine to the MODIS water- 
leaving radiance algorithm. 

• A modification to the SMA that does not require detailed aerosol models 
has been developed. This is important as the requirement for realistic 
aerosol models has been a weakness of the SMA. 

• We successfully acquired micro pulse lidar data in a Saharan dust outbreak 
during ACE-2 in the Canary Islands. 
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1. Atmospheric Correction Algorithm Development. 

a. Task Objectives: 

During CY 1997 there were seven objectives under this task. Task (i) below was considered 
to be the most critical. If the work planned under this task was successful, a module that enables 
the algorithm to distinguish between weakly- and strongly-absorbing aerosols would be included in 
the atmospheric correction algorithm. The seven tasks are listed below. 

(i) We will continue the study of the “spectral-matching” algorithm with the goal of having an 
algorithm ready for implementation by the end of CY 1997. 

(ii) We need to test the basic atmospheric correction algorithm with actual ocean color imagery. 
We will do this by looking at SeaWiFS and OCTS imagery as they become available. 

(in) We must implement our strategy for adding the cirrus cloud correction into the existing 
atmospheric correction algorithm. Specific issues include (1) the phase function to be used for 
the cirrus clouds, (2) the details of making two passes through the correction algorithm, and (3) 
preparation of the required look up tables (LUTs). 

(vi) The basic correction algorithm yields the product of the diffuse transmittance and the 
water-leaving reflectance. However, we have shown that the transmittance depends on the angular 
distribution of the reflectance only when the pigment concentration is very low and then only in the 
blue. We need to develop a method to include the effects of the subsurface BRDF for low-pigment 
waters in the blue. 

(v) We will initiate a study to determine the efficacy of the present atmospheric correction 
algorithm on removal for the aerosol effect from the measurement of the fluorescence line height 
(MOD 20). 

(vi) We will examine methods for efficiently including earth- curvature effects into the atmo- 
spheric correction algorithm. This will most likely be a modification of the look-up tables for the 
top-of-the-atmosphere contribution from Rayleigh scattering. 

(vii) We will examine the necessity of implementing our out-of-band correction to MODIS. 
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b. Work Accomplished: 

(i) We considered this task to be our most important atmospheric correction activity of 1997, 
and as such the major part of our effort on atmospheric correction was focussed on it. During this 
CY, we have further tested our spectral- matching algorithm (SMA) that, although very slow, is ca- 
pable of distinguishing between weakly- and strongly-absorbing aerosols. It is based on combining 
a model of the atmosphere with a water-leaving radiance model for the ocean, and effecting a vari- 
ation of the relevant parameters until a satisfactory fit to the MODIS top-of- atmosphere radiance 
is achieved. In simulations it showed significant success in detecting the absorption properties of 
the aerosol, i.e., distinguishing between weakly- and strongly-absorbing aerosols. We found ways to 
significantly increase the speed of the algorithm, and to enable it to operate using the same set of 
lookup tables that the basic algorithm uses. This enabled us to incorporate it in the basic correction 
algorithm, to be called each N x N pixels (where N ~ 10 — 100). In our present implementation the 
SMA simply returns a flag indicating the possible presence of absorbing aerosols. In subsequent im- 
plementations, we will return the oceanic and atmospheric properties. Unfortunately, at present we 
are not confident that models of absorbing aerosols, e.g., wind-blown dust, are sufficiently accurate 
to be effective candidate aerosol models, and this is impeding the full implementation. 

As with our original algorithm, the SMA requires candidate aerosol models, so it is important 
to understand the extent to which the candidate aerosol models need to be realistic. To try to 
understand this, we have revised the SMA by replacing the Shettle and Fenn aerosol models with 
Junge power-law models with a size-independent refractive index. This is an enormous simplifi- 
cation. We also replaced our “brute-force” search with a nonlinear optimization with which the 
aerosol size distribution, aerosol refractive index, and ocean properties are determined. Simulations 
show that excellent retrievals of the ocean’s pigment concentration are obtained; however, the price 
is a loss of accuracy in the estimate of the aerosol optical depth. A detailed report describing these 
results is presented as an appendix to this report. 

We have also expended much effort toward increasing the speed of the present MODIS algo- 
rithm. Significant gains have been achieved by not computing everything at each pixel. Rather, 
the LUT values that axe found for a given pixel are retained for 5-10 pixels. Although the is a small 
loss in accuracy, the increase in speed justifies it. This speed increase is an important, because the 
spectral-matching algorithm is very slow. 

(ii) Some imagery has been acquired from SeaWiFS. We have prepared a set of SeaWiFS- 
specific LUTs that will enable processing SeaWiFS imagery using the present MODIS code. This 
will allow a preliminary validation of the code. 
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(iii) None. In the light of the success of our spectral-matching algorithm, we may have to make 
significant modifications in our original strategy. This task has been put on hold to free resources 
for examination of task (i). The issues will be addressed during CY 1998 with the goal of having 
a complete implementation strategy ready during CY 1998. 

(iv) No work was carried out on this task. 

(v) No work was carried out on this task. 

(vi) No work was carried out on this task. 

(vii) The specifics of incorporating the out-of-band corrections in the MODIS algorithm have 
been worked out. 

c. Data/ Analysis/Interpretation: See item b above. 

d. Anticipated Future Actions: 

(i) We will continue work on the spectral-matching algorithm. Of particular interest is to 
implement the SMA optimization described in the Appendix. We will implement this method and 
test its efficacy with SeaWiFS imagery. We will examine the necessity of including bidirectional 
effects into the ocean model portion of the SMA. 

(ii) We shall continue testing the algorithm with SeaWiFS imagery through the first 7-9 months 
of CY 98, at which time we will switch to actual MODIS imagery. In particular, as described in 
(i), we want to test the SMA with real ocean color data. 

(iii) The cirrus cloud issue in the presence of our spectral-matching method needs to be 
explored. We will resolve the spectral-matching questions first, then devise a strategy to implement 
the cirrus correction. 

(iv) An ocean BRDF model will be developed in conjunction with the SMA (i), and validation 
(see 5 below). 

(v) We will perform a basic test of the efficacy of the correction algorithm for retrieving the 
fluorescence line height. 

(vi) None. 

(vii) None, until we are provided with the final MODIS spectral response functions. 
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f. Publications: 

H.R. Gordon, T. Zhang, F. He, and K. Ding, Effects of stratospheric aerosols and thin cirrus clouds 
on atmospheric correction of ocean color imagery: Simulations, Applied Optics , 36, 682-697 (1997). 

H.R. Gordon, Atmospheric Correction of Ocean Color Imagery in the Earth Observing System Era, 
Jour. Geophys. Res. y 102D, 17081-17106 (1997). 

H. Yang and H.E. Gordon, Remote sensing of ocean color: Assessment of the water-leaving radiance 
bidirectional effects on the atmospheric diffuse transmittance, Applied Optics , 36, 7887-7897 (1997). 

H.R. Gordon, T. Du, and T. Zhang, Remote sensing ocean color and aerosol properties: resolving 
the issue of aerosol absorption, Applied Optics, 36, 8670-8684 (1997). 
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2. Whitecap Correction Algorithm (with K.J. Voss). 

We believed that the basic objectives of this task had been realized, work was suspended until 
the validation phase, except insofar as the radiometer would be operated at sea when sufficient 
number of personnel were available. Our goal was to retain experience in operating and maintaining 
the instrumentation in preparation for the validation phase of the contract. However, by virtue of 
the surprisingly low reflectance increase due to whitecaps that we measured in the Tropical Pacific, 
we decided in July to reanalyze the data. This is a unique data set (our most important), as it 
was acquired in the trade winds with moderately high winds (8-12 m/s) and practically unlimited 
fetch and duration. 

a. Task Objectives: 

Operate the radiometer at sea to maintain experience in preparation for the validation phase 
and reanalyze the Tropical Pacific data. 

b. Work Accomplished: 

An undergraduate was hired part time and has learned to analyze the data. He is developing 
a differential filter to distinguish background from whitecaps, automatically. This is very tedious 
as it requires comparing the radiometer with the video monitor on a sample-to-sample basis. 

c. Data/Analysis/Interpretation 

At this time we have reduced the calibration data for the February cruise, but not the whitecap 
data. 

d. Anticipated Future Actions: 

First, we will reanalyze the Tropical Pacific whitecap data, then the data acquired during 
the February cruise. Second, we will operate the radiometer on the SeaWiFS initialization cruise 
in February 1998. This cruise will serve as a test bed for our initialisation cruise for MODIS in 
October. 

e. Problems/ Corrective Actions: None 

f. Publications: 

K.D. Moore, K.J. Voss, and H.R. Gordon, Spectral reflectance of whitecaps: Instrumentation, 
calibration, and performance in coastal waters, Jour. Atmos. Ocean. Tech. (In Press). 
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3. In-water Radiance Distribution (with K.J. Voss). 

a. Task Objectives: 

The main objective in this task is to obtain upwelling radiance distribution data at sea for 
a variety of solar zenith angles to understand how the water-leaving radiance varies with viewing 
angle and sun angle. 

b. Work accomplished: 

Data was acquired during the July 1997 cruise with this instrument. We have reduced this 
data to radiometric quantities and are now working with the data. During this cruise we had 
an opportunity to take an extensive data set during one afternoon (July 23), in which data was 
acquired with sun angles varying from approximately 40 degrees to 85 degrees. This will give us 

a good data set to evaluate the effect of sun angle on the upwelling radiance distribution at the 
MOBY site. 

c. Data/ Analysis/Interpretation: 


We have been concentrating our data analysis efforts on the extensive data set taken on July 
23. As an example, we show a polar-coordinate image of the upwelling radiance distribution at 
560 nm and a solar zenith angle of 66 degrees in the figure on the following page. Radiance 
increases from the red colors to the blue colors (blue highest). The colors are arranged linearly. 
This figure illustrates that the radiance distribution is neither uniform nor symmetric near the 
surface (the measurement depth was 2 meters). We are concentrating with our analysis on several 
cuts” through this data at 4 azimuthal angles. We will be comparing these images with results 
from a Monte Carlo model to test how well we can model these factors. In addition we will be 
comparing our measurements with the f /Q model of Morel (which is being used in other models 
in the community). 

d. Anticipated future actions: 

We will be continuing our data analysis. In addition we are calibrating the instrument and 
will be participating in another cruise with Dennis Clark in Hawaii during January and February 
of 1998, as well as the MODIS initalization cruise in October. 

e. Problems/Corrective actions: None. 

f. Publications: None. 
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4. Residual Instrument Polarization. 

The basic question here is: if the MODIS responds to the state of polarization state of the 
incident radiance, given the polarization- sensitivity characteristics of the sensor, how much will this 
degrade the performance of the algorithm for atmospheric correction, and how can we correct for 
these effects? 

a. Task Objectives: 

Add a module to perform the correction for residual instrument polarization. 

b. Work Accomplished: 

A module was added to perform the correction for residual instrument polarization. 

c. Data/Analysis/Interpretation: None. 

d. Anticipated Future Actions: 

Although this task is now basically complete. All that remains is incorporating the SBRS/MCST 
polarization-sensitivity data into the module. 

e. Problems/ Corrective Actions: None 

f. Publications: 

H.R. Gordon, T. Du, and T. Zhang, Atmospheric Correction of Ocean Color Sensors: Analysis of 
the Effects of Residual Instrument Polarization Sensitivity, Applied Optics , 36, 6938-6948 (1997). 
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5. Pre-launch/Post-launch Atmospheric Correction Validation (with K.J. Voss). 

a. Task Objectives: 

The long-term objectives of this task are four-fold: 

(i) First, we need to study aerosol optical properties over the ocean in order to verify the 
applicability of the aerosol models used in the atmospheric correction algorithm. Effecting this 
requires obtaining long-term time series in typical maritime environments. This will be achieved 
using a CIMEL sun/sky radiometer that can be operated in a remote environment and send data 
back to the laboratory via a satellite link. These are similar the radiometers used by B. Holben in 
the AERONET Network. 


(ii) Second, we must be able to measure the aerosol optical properties from a ship during 
the initialization/calibration/validation cruises. The CIMEL-type instrumentation cannot be used 
(due to the motion of the ship) for this purpose. The required instrumentation consists of an 
all-sky camera (which can measure the entire sky radiance, with the exception of the solar aureole 
region) from a moving ship, an aureole camera (specifically designed for ship use) and a hand- 
held sun photometer. We had a suitable sky camera and sun photometer but had to construct 
an aureole camera. Our objective for this calendar year was to make measurements at sea with 
this instrumentation, both to collect a varied data set and to test the instrumentation and data 
reduction procedures. 

In the case of strongly-absorbing aerosols, we have shown that knowledge of the aerosol vertical 
structure is critical. Thus, we need to be able to measure the vertical distribution of aerosols during 
validation exercises as well as build a climatology of the vertical distribution of absorbing aerosols. 
This will be accomplished with ship-borne LIDAR. We have procured a LIDAR system and modified 
it for ship operations. 

(iii) The third objective is to determine how accurately the radiance at the top of the atmo- 
sphere can be determined based on measurements of sky radiance and aerosol optical thickness at 
the sea surface. This requires a critical examination of the effect of radiative transfer on “vicarious” 
calibration exercises. 

(iv) The forth objective is to utilize data from other sensors that have achieved orbit (OCTS, 
POLDER, MSX), or are expected to achieve orbit (SeaWiFS) prior to the launch of MODIS, to 
validate and fine-tune the correction algorithm. 
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b. Work Accomplished: 

(0 The CIMEL was reinstalled during July in the Dry Tortugas and worked well for 2 months 
before failing in late September. We went down to remove it and send it back to Goddard for 
recalibration and found that it had fallen over (some bolts had been removed or not installed 
properly). The instrument was recalibrated and repaired at Goddard and we have reinstalled it 
in the Dry Tortugas (as of mid December). There were some problems during the recalibration, a 
filter was replaced and the tracking repaired, so the recalibration took longer than expected. We 
are investigating getting another CIMEL from the SeaWiFS instrument pool in the future. This 
will be used to swap out with our instrument during recalibration in the future, reducing downtime. 
In addition we are working on installing the Demonstrat program, which is used to view the data 
products, on one of our computers to allow easier access to the data. 

(ii) Aureole and sky camera data acquired during the July cruise was reduced during this 
period. These data were analyzed, specifically for several locations while the cruise went near the 
volcanic plume to look at the retrieved size distribution of particulates in the plume. 

The micro-pulse lidar system (MPL) was deployed in Tenerife, Canary Islands during June 
and July of 1997 in support of the Aerosol Characterization Experiment 2 (ACE-2). The primary 
goal of ACE-2 was to assess the radiative forcing due to anthropogenic aerosols over the North 
Atlantic. In addition, it was also possible to analyze the radiative impact of desert dust aerosols 
from North Africa. The primary purpose of operating the MPL during ACE-2 was to acquire 
information about the vertical structure of the dust layers. The second purpose for the ACE-2 
MPL deployment was to perform calibration procedures that are not possible in Miami and to test 
the MPL and planned lidar analysis algorithms for future MPL operations. Prior to the official start 
of ACE-2, an additional duty was assigned to the MPL operator by ACE-2 Headquarters. This 
consisted of faxing morning “real-time” lidar data, with initial estimates of the vertical location of 
any aerosol layers, to ACE-2 Headquarters for the morning preflight briefing of the ONR CERPAS 
Aircraft Pelican crew. The morning lidar data was used by the Pelican crew to aid in planning 
their flight paths. 

The MPL was shipped to Tenerife in May 1997 and was installed at Tenerife’s Observatorio de 
Izana (IZO) which is located on a mountain ridge near the center of the island at an altitude of 2367 
m above sea level. The location and altitude of IZO made it possible to perform both calibration 
procedures and study the vertical structure of dust layers moving west off the coast of Africa. The 
MPL was installed and began operational testing on June 14, 1997. Regular ACE-2 operation 
began on June 18, 1997. Regular MPL ACE-2 operation consisted of vertical, horizontal, and slant 
path (lidar tilted to a zenith angle of 60 degrees) measurements at specified times of each day. Slant 
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path measurements were made during midday to avoid direct sunlight from entering the detector 
and causing possible damage to the instrument. Vertical and horizontal measurements were made 
to study the vertical and horizontal structure of any aerosol layers and to perform calibration 
procedures. In addition, vertical measurements from 1000 GMT were used for the Pelican preflight 
plaining sessions. Regular MPL ACE-2 operation continued until shutdown on June 25, 1997. The 
MPL was then shipped back to Miami. 

(iii) The theoretical aspects of this work have been completed. The next phase is to use surface 
measurements to predict top-of-atmosphere radiance. 

(iv) We have prepared a duplicate version of the MODIS algorithm code to use the SeaWiFS 
spectral bands. This will be used to test the MODIS code with SeaWiFS data. 

c. Data/ Analysis/Interpretation: 

(ii) Two instrument problems occurred with the MPL during ACE-2. The first problem was 
due to a poor cable connection from the lidar receiver to the MPL scalar unit. Normal nighttime 
temperatures at IZO during June were around 10°C. Contraction of the BNC cable connector due 
to the cold weather resulted in the loss of nighttime data from June 18 to June 27, 1997. This 
problem was fixed by changing cables on June 28, 1997. The second MPL problem to occur during 
ACE-2 was attributed to the MPL laser. The laser doubling crystal burned midway through ACE-2 
due to a poor ground connection to the laser temperature controller. This resulted in distortion of 
the outgoing laser pulses and a gradually decreasing signal- to- noise ratio throughout the month of 
July. It was not possible to perform normal MPL calibration procedures for the ACE-2 MPL data 
because of the laser problem. A modified calibration procedure was developed which normalized 
the lidar signals to a Rayleigh only lidar signed using MPL data from clean nighttime periods. This 
modified calibration procedure was used to correct all the ACE-2 MPL data. 

The initial results using the modified calibration procedure are very promising. In particular, 
the day of July 17, 1997 has been analyzed to study the validity of the calibrations. There were 
three significant dust events during ACE-2, the second occurred during the period July 16 - July 18, 
1997. The figure on the following page shows a time series plot of attenuated backscatter measured 
with the MPL during this dust event. The rainge-corrected lidar signal is defined as: 

S(z) = CE[f3 R (z) + 0a{*)] exp[-2[r fl (z) + r x (z)]], 

where z is the range (km), C is the instrument constant, E is the output energy, f3(z) is the 
backscatter coeff. (1/km), and r(z) is the opticad depth. The subscripts R and A denote Rayleigh 
and Aerosol respectively. Attenuated backscatter is defined as S(z)/CE. It is the totail backscatter 
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Figure 1: Attenuated Backscatter Coeff. (July 16 - July 18) 
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coefficient times the transmission squared. This quantity is not dependent on an inversion algorithm 
and is the initial data product from the MPL. The figure shows the beginning of the dust event 
(approximately 2100 GMT on July 16, 1997) as well as the vertical structure of the layer until its 
disappearance during the evening of July 18, 1997. Gaps in the figure are due to periods where 
the MPL was performing horizontal or slant path measurements. Horizontal and slant path ACE-2 
data have not yet been analyzed. 

On July 17, 1997 at 1830 GMT, simultaneous measurements of the aerosol optical depth 
(AOD) were made at IZO by a NASA AERONET Cimel sunphotometer 1 , as well as on top of Teide 
(Tenerife’s volcano, altitude 3500 m) with a YES shadowband 2 , and by the Pelican’s suntracking 
sunphotometer 3 as the aircraft spiraled down above IZO from an altitude of approximately 3800 
m. Recent comparisons of the MPL AOD vertical profile with the AOD profile of from Pelican 
and the AOD from Teide show very good agreement (within 0.03 AOD units at all altitudes from 
2367 m to 3800 m). The AOD from the IZO Cimel was used as an input constraint in the MPL 
inversion algorithm and thus was not part of the lidar comparison, however, comparison of the 
IZO Cimel AOD with the Pelican AOD at 2367 m agree within about 0.01 AOD units. The good 
agreement between the MPL AOD profiles and those of independent instruments, as well as analysis 
of corrected clean nighttime lidar data (inversions of corrected data do show clean, insignificant 
aerosol amounts) show that the modified correction procedure has been a success. 

d. Anticipated Future Actions: 

(i) The CIMEL is reinstalled in the Dry Tortugas, and we Me looking forward to collecting 
coincident data from the CIMEL and SeaWiFS satellite. In particular we would like to see data 
like this after a winter cold front, when the atmosphere is typically clearer than normal. 

(ii) The aureole and sky camera are being recalibrated in anticipation of another cruise in 
January- February. We will also continue with the data reduction for these instruments. 

Future use of the ACE-2 MPL data is currently being planned. The data will be made directly 
available to all ACE-2 participants as outlined in the ACE-2 guidelines. The data will also be 
the focus of at least two planned papers. The first will describe the MPL ACE-2 system and the 
problems and resulting calibration procedures used to correct the data. The second paper will 
involve other ACE-2 participants and will focus on closure analysis of the type described above for 
the day of July 17, 1997. The comparison of lidar derived aerosol optical properties, such as the 

1 B. Holben, A. Smirnov (personal communication) 

2 P. Formenti, M. Andreae (personal communication) 

3 B. Schmid, P. Russell, J. Livingston, P. Durkee (personal communication) 
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extinction coefficient and the AOD, with those of other independent instruments in the near-range 
lidar field (the lidar range most effected by lidar overlap and afterpulse problems, 0-5 km) is a 
novel and important comparison in the field of lidar research. The initial agreement between the 
MPL AOD and the other instruments (within 0.03 AOD units) demonstrates the possibility of 
using MPL data to acquire accurate aerosol vertical profiles. In addition to the two planned ACE-2 
papers, the ACE-2 MPL data will be a primary portion of the dissertation of graduate student E.J. 
Welt on. 

The MPL has been repaired and returned to Miami. We will deploy it during the January - 
February SeaWiFS initialization cruise and the October-November MODIS initialization cruise. 

e. Problems/corrective actions: None. 

f. Publications: 

T. Zhang and H.R. Gordon, Columnar aerosol properties over oceans by combining surface and 
aircraft measurements: sensitivity analysis, Applied Optics , 36, 2650-2662 (1997). 

T. Zhang and H.R. Gordon, Retrieval of elements of the columnar aerosol scattering phase matrix 
from sky radiance over the ocean: simulations, Applied Optics , 36, 7948-7959 (1997). 

D.K. Clark, H.R. Gordon, K. J. Voss, Y. Ge, W. Broenkow, and C. Trees, Validation of Atmospheric 
Correction over the Oceans, Jour. Geophys. Res., 102D 17209-17217 (1997). 

K.J. Voss and Y. Liu, Polarized radiance distribution measurements of skylight: Part 1, system 
description and characterization, Applied Optics 36, 6083-6094 (1997). 

Y. Liu and K.J. Voss, Polarized radiance distribution measurements of skylight: Part 2, experiment 
and data, Applied Optics 36, 8753-8764 (1997). 

H. Yang and H.R. Gordon, Retrieval of the Columnar Aerosol Phase Function and Single Scattering 
Albedo from Sky Radiance over Land: Simulations, Applied Optics (Accepted). 

H.R. Gordon, In-orbit calibration strategy of ocean color sensors, Remote Sensing of Environment 
(Accepted). 
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6. Detached Coccolith Algorithm and Post Launch Studies (W.M. Balch). 

a. Task Objectives: 

Our MODIS work during involves understanding all aspects of the influence of suspended cal- 
cium carbonate particles on inherent optical properties in the sea. Work during this reporting 
period focused on several areas: processing cell count, atomic absorption, particulate organic car- 
bon, and scattering data from MODIS pre-launch cruises in the Gulf of Maine, preparation and 
submission of three MODIS-related manuscripts, another MODIS pre-launch cruise to the Gulf of 
Maine to examine in situ properties of cal cite particles, and analysis calcite optical properties from 
the Arabian Sea. 

The algorithm for retrieval of the detached coccolith concentration from the coccolithophorid, 
E. huxleyi is described in detail in our ATBD. The key is quantification of the back scattering coef- 
ficient of the detached coccoliths. Our earlier studies focused on laboratory cultures to understand 
factors affecting the calcite-specific backscattering coefficient. A thorough understanding of the 
relationship between calcite abundance and light scattering, in situ, will provide the basis for a 
generic suspended calcite algorithm. As with algorithms for chlorophyll, and primary productivity, 
the natural variance between growth related parameters and optical properties needs to be under- 
stood before the accuracy of the algorithm can be determined. To this end, the objectives of our 
coccolith studies during this reporting period have been: 

(1) acquire optical field data on the distribution and abundance of coc- 
colithophores in the Gulf of Maine, 

(2) process data from recent cruises, and 

(3) publish earlier results from our remote sensing work done with the 
CalCOFI data set, flow cytometry experiments, optical experiments 
with coccoliths, and Gulf of Maine pre-launch cruises. 

For perspective on the directions of our work, we provide an overview of our previous activities: 

Jan-June 1995: Research focus - chemostat cultures (in which algal 
growth rate was precisely controlled) and we examined how the optical 
properties of these calcifying algae changed as a function of growth. 
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July-Dee 1995: Research focus - shipboard measurements of suspended 
calcite and estimates of optical backscattering as validation of the lab- 
oratory measurements. We participated on two month-long cruises to 
the Arabian sea, measuring coccolithophore abundance, production, and 
optical properties. 

Jan-June 1996: Research focus - field calcite distributions, during two 
Gulf of Maine cruises, one in March and one in June. 

July-Dee 1996: Research focus - participated on another cruise to the 
Gulf of Maine and processed samples from the Gulf of Maine. 

Jan-June 1997: Research focus - continued processing samples from our 
previous cruises, upgraded our laser light scattering photometer used in 
all of the calcite scattering measurements, performed another pre-launch 
cruise on calcite particle optics in the Gulf of Maine, and analyzed our 
results from the MODIS-funded flow cytometer work. 

b. Work Accomplished: 

1) Processing of the suspended calcite samples from the June ’97 cruise. We processed the 
CHN samples from this last trip, which will be important for understanding the backscattering 
arising from suspended organic matter. 

2) Final analysis and write-up of our flow cytometer experiments. In October ’97, this 
manuscript was submitted for the MODIS special issue to appear in J. Geophys. Res. next year. 

3) Continued microscope cell/coccolith counts for samples from the Gulf of Maine, We are 
processing count data from ’96, and working on enumeration of ’97 samples. 

4) A paper on coccolith optics with Ken Voss and Katherine Kilpatrick, has been reviewed, 
revised, and sent back to Limnology and Oceanography. It should appear in print in 6-9 months. 
The paper deals with a laboratory comparison of volume scattering instruments, using the coccol- 
ithophore, Emiliania huxleyi, and how to model the optical properties. 

5) We performed a pre-launch MODIS cruise in the Gulf of Maine aboard the R/V Delaware 
II, and as with our previous Albatross cruise, weather conditions were extremely bad (up to 50kt 
winds and 20 foot seas). We ran our underway system for the entire trip (~ 1700 nautical miles) 
without a problem. On this cruise, we provided sea truth radiance and irradiance data for the 
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SeaWiFS (with imagery coordinated with Dr. Gene Feldman). We visited 60 stations, 28 of which 
were full optical stations (with measurements of total backscattering, acid-labile back scattering, 
calcite concentrations, coccolithophore and coccolith counts). We had the Yentsch/Phinney group 
on board measuring spectral absorption, in situ absorption/attenuation (AC-9), upwelling and 
downwelling irradiance (S ATLANTIC TSRB bio-optical sampler), and dissolved organic matter 
absorption. We had 5 days with clear conditions, and radiometers in the water for the SeaWiFS 
overpass. 

6) Suspended calcite samples from the Gulf of Maine have been run in the graphite furnace 
atomic absorption spectrometer at the University of Maine. We now only have samples from our 
most recent cruise in November ’97. 

7) All cell and coccolith counts from the Arabian Sea are now being analyzed for publication. 
The first sections of coccolith concentration are now complete, and show some fascinating trends. 
Data from JGOFS Process Cruise 6 are now entered into spreadsheets. In terms of current mi- 
croscopy, we now are focusing on the several hundred cell count samples from our 1997 Gulf of 
Maine cruises. This will take several months to complete. 

8) Calcification data from the March, June and November 1996 Gulf of Maine cruises have been 
analyzed and the first draft of a manuscript is now complete. The manuscript will be submitted 
for publication following final changes. 

c. Data/ Analysis/Interpretation: 

Single cell experiments : 

The results of our calcite flow cytometry work have been submitted for publication in the 
MODIS special issue of JGR-Oceans. 

Cruise results : 

As with our June ’97 cruise, calcite-dependent backscattering was quite high in the Gulf of 
Maine during November, 1997. This cruise was the first late Fall cruise in the Gulf of Maine, where 
calcium carbonate optical properties were quantified. Calcite scattering commonly accounted for 
10-20% of total backscattering. Interestingly, even with the low light conditions, coccolithophores 
were still remarkably abundant in the Gulf of Maine. We formerly interpreted strong coccolith 
influence on particle optics as a H summer-only”issue but it now can be clearly considered a year- 
round phenomenon. 
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Given the future plan for more Gulf of Maine cruises, as well as our SIMBIOS activities, we are 
writing new software to merge the various data sets into one file that can be submitted to NASA’s 
Sea Bass data archive. The first version of the software (which performs vicarious calibrations, offset 
corrections, plus calculating “derived quantities”^ now complete. This software will considerably 
reduce the time currently required to quality control the numbers and produce hydrographic plots 
of the data. 

The March, June and November cruises to the Gulf of Maine were processed to show the aerial 
distribution and depth profiles of calcification, primary productivity, and calcite concentrations. 
These results represent the largest assemblage of their kind in the Gulf of Maine, and are featured 
in the manuscript currently in draft form. 

d. Anticipated Future Actions: 

Work in the next year will address severed areas: 

1) Processing of the suspended calcite samples from the November ’97 cruise. 

2) Any necessary revision of our flow cytometer manuscript. 

3) Continued microscope cell/coccolith counts for samples from the Gulf of Maine. 

4) We will go to sea in June ’98 on yet another MODIS Gulf of Maine cruise. 

e. Problems/Corrective Actions: None 
f* Publications: 

Voss, K., W. M. Balch, and K. A. Kilpatrick. Scattering and attenuation properties of Emiliania 
huxleyi cells and their detached coccoliths. Limnol. Oceanogr. (In press). 

Balch, W. M. and B. Bowler. Sea surface temperature gradients, baroclinicity, and vegetation 
gradients in the sea. J. Plank. Res . (In press). 

Balch, William M., David T. Drapeau, Terry L. Cucci, and Robert D. Vaillancourt, Katherine 
A. Kilpatrick, Jennifer J. Fritz. Optical Backscattering by Calcifying Algae. J. Geophys. Res. 
(Submitted). 
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Campbell, J., D. Antoine, R. Armstrong, W. Balch, R. Barber, and others. Comparison of algo- 
rithms for estimating primary productivity from surface chlorophyll temperature and irradiance. 
Global Biogeochem . Cycles (In preparation). 

Graziano, L., W. Balch, D. Drapeau, B. Bowler, and S. Dunford. Organic and inorganic carbon 
production in the Gulf of Maine. Coni. Shelf Res, (In preparation). 

Fritz, J. and W. Balch. Verification of carbon- 14 calcification rates of the coccolithophore Emiliania 
huxleyi using scanning electron microscopy. (In preparation). 
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7. Other Activities. 


The PI participated in the SIMBIOS meeting in Maryland in August, the MODIS Science Team 
Meeting in October, and a SeaWiFS data acceptance review in November. K. Voss participated in 
an INDOEX planning meeting in September, and a SIMBIOS remote-sensing-reflectance workshop 
in Norfolk in December. 

Two papers were presented at the AGU Fall Meeting in San Francisco. They are provided 
below. 

J. Ritter, K. Voss, H. Gordon, D. Clark, and L. Koval, Solar Aureole Measurements of Volcanic 
Aerosol in the Kilauea Plumb, EOS Transactions , 78-46, A41A-27. 

E. J. welton and K.J. Voss, Preliminary Lidar Results from the Aerosol Characterization Experiment— 
2 (ACE-2), EOS Transactions , 78-46, A41A-13. 
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8. 1997 Publications. 

H.R. Gordon, T. Zhang, F. He, and K. Ding, Effects of stratospheric aerosols and thin cirrus clouds 
on atmospheric correction of ocean color imagery: Simulations, Applied Optics, 36, 682-697 (1997). 

H.R. Gordon, Atmospheric Correction of Ocean Color Imagery in the Earth Observing System Era, 
Jour. Geophys. Res., 102D, 17081-17106 (1997). 

H. Yang and H.R. Gordon, Remote sensing of ocean color: Assessment of the water-leaving radiance 
bidirectional effects on the atmospheric diffuse transmittance, Applied Optics, 36, 7887-7897 (1997). 

H.R. Gordon, T. Du, and T. Zhang, Remote sensing ocean color and aerosol properties: resolving 
the issue of aerosol absorption, Applied Optics, 36, 8670-8684 (1997). 

H.R. Gordon, T. Du, and T. Zhang, Atmospheric Correction of Ocean Color Sensors: Analysis of 
the Effects of Residual Instrument Polarization Sensitivity, Applied Optics, 36, 6938-6948 (1997). 

T. Zhang and H.R. Gordon, Columnar aerosol properties over oceans by combining surface and 
aircraft measurements: sensitivity analysis, Applied Optics, 36, 2650-2662 (1997). 

T. Zhang and H.R. Gordon, Retrieval of elements of the columnar aerosol scattering phase matrix 
from sky radiance over the ocean: simulations, Applied Optics, 36, 7948-7959 (1997). 

D.K. Clark, H.R. Gordon, K.J. Voss, Y. Ge, W. Broenkow, and C. Trees, Validation of Atmospheric 
Correction over the Oceans, Jour. Geophys. Res., 102D 17209-17217(1997). 

K.J. Voss and Y. Liu, Polarized radiance distribution measurements of skylight: Part 1, system 
description and characterization, Applied Optics 36, 6083-6094 (1997). 

Y. Liu and K.J. Voss, Polarized radiance distribution measurements of skylight: Part 2, experiment 
and data, Applied Optics 36, 8753-8764 (1997). 

K.D. Moore, K.J. Voss, and H.R. Gordon, Spectral reflectance of whitecaps: Instrumentation, 
calibration, and performance in coastal waters, Jour. Atmos. Ocean. Tech. (In Press). 

H.R. Gordon, In-orbit calibration strategy of ocean color sensors, Remote Sensing of Environment 
(In Press). 

Balch, W. M. and B. Bowler. Sea surface temperature gradients, baroclinicity, and vegetation 
gradients in the sea. J. Plank. Res. (In Press). 
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H. Yang and H.R. Gordon, Retrieval of the Columnar Aerosol Phase Function and Single Scattering 
Albedo from Sky Radiance over Land: Simulations, Applied Optics (Accepted). 

Voss, K., W. M. Balch, and K. A. Kilpatrick. Scattering and attenuation properties of Emiliania 
huxleyi cells and their detached coccoliths. Limnol. Oceanogr. (In press). 

Balch, W. M. and B. Bowler. Sea surface temperature gradients, baroclinicity, find vegetation 
gradients in the sea. In press, /. Plank. Res. (In press). 

Balch, William M., David T. Drapeau, Terry L. Cucci, and Robert D. Vaillancourt, Katherine 
A. Kilpatrick, Jennifer J. Fritz. Optical Backscattering by Calcifying Algae. J. Geophys . Res. 
(Submitted). 
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Appendix 


Atmospheric correction of ocean color imagery: 

Use of the Junge power-law aerosol size distribution with 
variable refractive index to handle aerosol absorption 



Abstract 


When strongly- absorbing aerosols are present in the atmosphere, the usual two-step 
procedure of processing ocean color data — (1) atmospheric correction to provide the 
water-leaving reflectance {pw)^ followed by (2) relating p w to the water constituents — 
fails and simultaneous estimation of the ocean and aerosol optical properties is necessary. 
We explore the efficacy of using a very simple model of the aerosol — a Junge power-law 
size distribution consisting of homogeneous spheres with arbitrary refractive index — in 
a nonlinear optimization procedure for estimating the relevant oceanic and atmospheric 
parameters for Case 1 waters. Using simulated test data generated from more realis- 
tic aerosol size distributions (sums of log-normally distributed components with different 
compositions), we show that the ocean’s pigment concentration ( C ) can be retrieved with 
good accuracy in the presence of weakly- or strongly-absorbing aerosols. However, because 
of significant differences in the scattering phase functions for the test and power-law distri- 
butions, large error is possible in the estimate of the aerosol optical thickness. The positive 
result for C suggests that the detailed shape of the aerosol scattering phase function is not 
relevant to the atmospheric correction of ocean color sensors. The relevant parameters are 
the aerosol single scattering albedo and the spectral variation of the aerosol optical depth. 
We argue that the assumption of aerosol sphericity should not restrict the validity of the al- 
gorithm, and suggest an avenue for including colored aerosols, e.g., wind-blown dust, in the 
procedure. A significant advantage of the new approach is that realistic multi-component 
aerosol models are not required for the retrieval of C . 
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1. Introduction 


It is now generally recognized that measurement of the radiance reflected from the 
ocean atmosphere system can be used to estimate the concentration of phytoplankton pig- 
ments, the sum of the concentrations of chlorophyll a and phaeophytin a, in the ocean. 1-3 
(At present, this estimation is generally limited to Case 1 waters, 3 i.e., waters in which 
phytoplankton and their immediate detrital material constitute the variable component of 
the water’s optical properties.) Thus several sensors have been, or are soon to be, launched 
for the purpose of studying the spatial-temporal variation of phytoplankton globally, to 
understand the role of phytoplankton and their variability in global climate. The sea- 
viewing wide-field-of-view sensor 4 (SeaWiFS) was launched on August 1, 1997, and the 
moderate resolution imaging spectroradiometer 5 (MODIS) is scheduled for launch in the 
summer of 1998. In Case 1 waters, the component of the measured radiance that carries 
information about the phytoplankton concentration, the radiance backscattered out of the 
water and transmitted to the top of the atmosphere (the water-leaving radiance), is at 
most 10% in the blue and typically much smaller in the green. The rest of the signal is 
comprised of radiance reflected from the atmosphere and the sea surface. Thus, the water- 
leaving radiance must be extracted from the measured radiance — a process referred to 
as atmospheric correction. 

Typically, atmospheric correction algorithms assume that radiative transfer in the 
ocean and atmosphere can be decoupled, i.e., photons are never backscattered out of 
the ocean more than once. 6 Furthermore, determination of the phytoplankton pigment 
concentration is carried out in two steps: (1) the water-leaving radiance is retrieved as a 
product of the atmospheric correction; and (2) the pigment concentration is determined 
from the water-leaving radiance. The atmospheric correction algorithm has increased in 
complexity from the simple single-scattering algorithm 7 ’ 8 used for the coastal zone color 
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scanner (CZCS) — the proof-of-concept ocean color mission — to the complete multiple- 
scattering algorithm proposed for SeaWiFS. 9 ’ 10 The new algorithm is capable of providing 
the water-leaving radiance with the required accuracy for typical marine atmospheres; 
however, it fails when the atmosphere contains strongly- absorbing aerosols, e.g., in the 
presence of urban pollution or desert dust transported over the ocean by the winds. 

Recently, Gordon, Du, and Zhang 11 proposed an algorithm for Case 1 waters that is 
capable of atmospheric correction in the presence of both weakly- and strongly-absorbing 
aerosols. In this algorithm, the aerosol properties and the pigment concentration are de- 
termined simultaneously. This simultaneous determination is required because the aerosol 
absorption is manifest mostly in the blue, where the water-leaving radiance is usually 
largest, and the most variable. In effect, in this algorithm the ocean pigment concentra- 
tion and the aerosol model are systematically varied and the top-of-the-atmosphere (TOA) 
radiance is computed. The candidate aerosol models used in this algorithm were derived 
from the bimodal log-normal aerosol size distributions and refractive indices provided by 
Shettle and Fenn 12 as a function of relative humidity (RH). The water-leaving radiance 
was computed as a function of the pigment concentration (C) and a marine particulate 
scattering factor (6°) using the semi-analytic model of Gordon et al. 13 The mean of the 
ocean-atmosphere properties that yield the ten-best fits (in an RMS sense) of the com- 
puted to the measured radiances provides the derived parameters (pigment concentration, 
aerosol optical thickness, aerosol single scattering albedo, etc.). 

Zhao and Nakajima 14 proposed a similar algorithm for CZCS, i.e., the aerosol prop- 
erties and the pigment concentration determined simultaneously; however, because of the 
severe limitation of the number of spectral bands on CZCS, they employed models with 
Junge 15 power-law size distributions with a single index of refraction, and used ratios of 
water-leaving radiances rather than the water-leaving radiance itself. They showed good 
retrievals with simulated test data created using bimodal aerosol size distributions with 
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the same index of refraction as the Junge power-law distribution, and concluded that the 
power-law distribution was an adequate approximation to the now-favored bimodal aerosol 
distributions for the purposes of atmospheric correction. In contrast to bimodal log-normal 
size distributions, which require several parameters to describe (the modal diameters and 
standard deviations of both modes along with their relative concentrations), the simplicity 
of the Junge power-law distribution is attractive in that only one parameter is required. 
A further simplification is that the bimodal distributions require an index of refraction for 
each mode (typically different), while a single refractive index is usually assumed for the 
power-law distribution. 

In this paper, we explore replacing the Shettle and Fenn 1 Candidate aerosol bimodal 
models in the Gordon, Du, and Zhang 11 procedure with a Junge power-law model; however, 
in order to be able to deed with both weakly- and strongly-absorbing aerosols, we let the 
real and imaginary parts of the index of refraction vary as well as the size distribution 
parameter — a natural extension of the Zhao-Nakajima 14 algorithm. Thus, all of the 
properties of the aerosol are allowed to vary continuously, subject to the constraint of the 
given mathematical form for the aerosol size distribution. This allows use of nonlinear 
optimization techniques to effect a solution. The results show that excellent retrievals of 
the pigment concentration are obtained, weakly- and strongly-absorbing aerosols are easily 
identified, but the retrieved aerosol optical thickness can have large errors because of the 
significant differences in the scattering phase functions for bimodal log-normal and power- 
law distributions. Thus, we also answer the question, is a correct aerosol phase function 
really required for atmospheric correction of ocean color imagery? 

We begin with a description of the aerosol and ocean models used in the procedure and 
describe the nonlinear optimization. Next, we provide tests with simulated data that is free 
of radiometric calibration errors, and follow with an examination of the effect of realistic 
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radiometric calibration errors. Finally, we examine the performance of the algorithm as a 
function of the aerosol’s vertical distribution. 

2. The correction mechanics 

We begin by using a generalized reflectance p in place of radiance L. These are related 
by p = ttL/Fq cos do, where Fq is the extraterrestrial solar irradiance and do is the solar 
zenith angle. Then, after avoiding sun glitter and correcting for the presence of whitecaps, 
the remaining upwelling reflectance at the top of the atmosphere p t (A) consists of the 
following components: 10 the pure Rayleigh (molecular) scattering contribution p r (A) , the 
pure aerosol scattering contribution p a (A), the contribution due to the interaction effect 
between air molecules and aerosols p ra (A), and the desired water-leaving contribution 
*(A)/>w(A), i.e., 

/><( A) = p r ( A) + p a ( A) + p r a{ A) + t(A)p l „(A), (1) 

where <( A) is the diffuse transmittance of the atmosphere. As p r ( A) can be computed 
precisely from an estimate of the surface atmospheric pressure, 16 ’ 17 it can be subtracted 
from p t ( A) to form 

[pt(A) — p r (A)] = [p a ( A) +Pra(A)] + [f(A)p ti ,(A)j . (2) 

Following Gordon, Du, and Zhang, 11 we use an atmospheric model to compute p a ( A) + 
p ra ( A) and f(A), an ocean model to compute Pu>(A), and compare the results to the given 
value of p t ( A) — p r (A). 

A. The water and aerosol models 
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The water reflectance model is derived from Gordon et al., 13 and is identical to that 
employed by Gordon, Du, and Zhang. 11 It has two parameters C , the pigment concen- 
tration, and 6° a parameter related to particulate scattering. 3 The model assumes that 
the upwelling spectral radiance below the sea surface is totally diffuse, contrary to the 
observations of Morel and co-workers. 18-21 However, as Morel and Gentili 21 have shown, 
bidirectional effects can be easily introduced into the model and described as a function of 
C (and, if necessary, 6°). 


The aerosol model utilizes a Junge power-law size distribution: 


dN 

Id 
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Do <D < D \ , 
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= 0, 


D > D 2} 


where dN is the number of particles per unit volume with diameters between D and 
D + dD. In this paper, we choose D 0 = 0.06 /mi, D x = 0.20 /zm, and D 2 — 20 /zm. 
The size distribution is thus characterized by a single parameter, v. A refractive index 
m = m r — irrii is then combined with this size distribution to compute the aerosol optical 
properties — the phase function, optical thickness (r a (A)), and single scattering albedo 
(u;o(A)) — using Mie theory. As the synthetic test data is also computed from bimodal 
aerosol models using Mie theory (Subsection 3A), this procedure does not introduce any 
error, i.e., we assume all particles are spherical. The aerosol optical properties are then 
used to generate the reflectance p a (A) + />ra(A) as a function of the aerosol optical thickness 
for all viewing-sun geometries. A two-layer radiative transfer code (aerosols in the lower 
layer) in which the air-sea interface is modeled by a flat Fresnel-reflecting surface, is used 
in the computations. These reflectances are stored in the form of lookup tables (LUTs) 
for rapid access by the algorithm. In our LUTs, the values of the parameter v range from 
2 to 4.5 in steps of 0.5, and m* = 0, 0.001, 0.003, 0.01, 0.03, and 0.04, with m r = 1.333 
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and 1.50. Thus, we have a total of 72 distinct power-law aerosol models that can be used 
to interpolate for arbitrary values of z/, m r , and m;. We use linear interpolation for u and 
quadratic interpolation for m T and ra; [Formula 25.2.66 of Ref. 22]. 


In the optimization procedure, there are a total of six parameters: C and b° for the 
water, v, m r , mi, and r a (865), for the atmosphere. [The aerosol optical thickness at 865 
nm, r a (865), is the measure of the aerosol concentration.] However, we can reduce the 
number of parameters by noting that p w % 0 in the near infrared (NIR). [SeaWiFS has 
NIR bands at 765 and 865 nm, while MODIS has bands at 748 and 869 nm.] For these 
bands, the value of p a (A) + p r a(A) is known. The spectral variation of p a (A) + />ra(A) in the 
NIR has been shown to depend mostly on the aerosol size distribution and only slightly on 
the refractive index, 10 thus, these two bands can be used to determine the size distribution 
parameter v. This is effected by comparing the measured value of 


e'(A„A<) 


Pa(A a ) + /?ra(A$) 

Pa(A^) + Pra(A^) 


( 3 ) 


where X s and A* are the shorter and the longer of the two NIR bands available on the sensor, 
with that computed from the models as a function of v. Figure 1 shows an example of the 
relationship between v and £'(765,865) for the power-law distributions. The figure was 
prepared using the 12 values of the complex refractive index and four values of r a (865). It 
shows that a coarse estimate of v (maximum error ~ +0.3 to —0.1, standard error ~ +0.15 
to —0.1) can be obtained from the measured £'(765,865). Given the value of v determined 
from the i/-e f (765, 865) relationship for the given geometry, it is a simple matter to compute 
the aerosol optical thickness that provides the correct value of reflectance at A t for each 
m r — irrii. Thus, both v and r a (865) can be removed from the list of unknown parameters 
via the NIR measurements. In the case of SeaWiFS, this leaves six bands from which the 
four remaining parameters must be determined. It should be noted that using this method 
of reducing the number of parameters introduces some error. For example, if the actual 
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aerosol model is identical to one of the 72 candidate models, a close, but incorrect, value 
of v will usually be chosen by the algorithm. 

B. The nonlinear optimization. 

Before continuing the discussion we introduce some simplifying notation: 

PA ! TYli , 77Z r ) = Pa (A* 5 TTl { , 77i r ) 4“ pra (Ai ! 771 r ) 

(4) 

Pw = t(\i) p w (Ai ; (7, 6°), 

where p' A is the resulting reflectance of the aerosol and interaction term computed at 
A; for the Junge power-law size distribution with the parameters v and r a (865) (both 
determined as in Subsection 2A) using m* and m r , and p* w is the diffusely transmitted 
normalized water-leaving reflectance computed at A; for water parameters C and 6°. The 
assumption that p\ (A; and p ! w (A; ; (7, 6°) form the final reflectance field on the 

top of the atmosphere (TOA) independently leads to the system of nonlinear equations in 
four variables (mi, m r , (7, 6°): 

f Pa(*i ;wii,m r ) +p' w (Xi -,C,b°) = /? t (Ai) - p r (\i) 

p' A (A 2 ; rrij, m r ) + p w (A 2 \C,b°) = p t (A 2 ) - />r(A 2 ) 

< (5) 

, Pa{^n v + p' w (\ Nv ;C,b°) — pt(>^N v ) ~ pr(^N v )- 

We seek to solve this set of equations for mi, m r , C , and b° when N v , the number of 

visible bands used in the retrieval, is greater than the number of unknowns (4), and adopt 
a nonlinear optimization procedure in order to do so. We note that the set of variables 
(■mi,m r ,C, b°) is subject to the constraints: 

' 0 < mi(min) < rrii < mj(max) 

0 < m r (min) < m r < m r (maz) 

< ( 6 ) 

0 < C(min) < C < C(max ) 

k 0 < b°(min) < 6° < b°(max). 
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As in Gordon, Du, and Zhang, 11 we assume that C(min ) = 0.05 mg/m 3 , C(max) — 1.50 
mg/m 3 , b°(min ) = 0.12 m -1 , and b°(max) = 0.45 m -1 . For the atmosphere, we take 
mi{rnin ) = 0, rrii(max ) = 0.04, m r (min ) = 1.33, and ra r (mai) — 1.50, as realistic ranges 
for the aerosol parameters. We choose as the the least-square objective function: 


£.2 / f, h Qx 1 V- / p\ (•'i (Ai ;C,6°) V , . 

S lsQ (m i ,m.,C.4) = p ? - rT )L( k rtV-rtt) )' () 


At the solution point (m*, m*, C*, b 0 *), Slsq should attain its minimal value subject to 
constraints Eq. (6). 


A number of versatile methods exists to find such a solution 23 although they might 
differ appreciably as to the convergence rate, the memory requirements, and the ability to 
converge to the solution. In particular, methods belonging to the so-called quasi-Newton 
class, 24 that take advantage of the information contained in the Hessian matrix of a func- 
tion (i.e., the square matrix of the second partial derivatives of the function evaluated at a 
point) in calculating the step size and the direction of search, appear to be adequate for our 
purposes. Powell 25 proves the convergence to a minimum for the Davidon-Fletcher-Powell 
(DFP) variable metric class method [which also holds for Broyden-Fletcher-Goldfarb- 
Shanno (BFGS) algorithm] on a function with properties outlined above [Eqs. (5) and 
(6)]. However, to the best of our knowledge, the quasi-Newton class methods do not guar- 
antee, in general, a convergence to a minimum. In this study we used DZXMWD routine 
from the commonly available IMSL library, 26 which is based on the Harwell library routine 
VA10A and incorporates a quasi-Newton method. To assure the convergence to a mini- 
mum in our non-linear optimization procedure a set of starting search points is taken, some 
of which are discarded after performing a few iterations based on their relative values. We 
observed that the resulting solution was the same (within a pre-set tolerance) whatever 
the number of starting search points were chosen, meaning that convergence to a global 
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minimum was probably achieved, so we believe it is sufficient to use just one starting point 
in solving the system [Eq. (5)] and the result is the solution we want. 

3. Algorithm performance with simulated data. 


A. Characteristics of the simulated p t ( A) data. 


As in earlier work, 9 ’ 11 we tested the performance of the algorithm with simulated data 
— pseudodata. The pseudodata were created by choosing specific values of C and b° for 
the computation of p w ( A), and by utilizing specific aerosol models for the computation of 
Pr Pa d - Pro. • The aerosol models we chose for the tests were those developed by Shettle and 
Fenn 12 at a relative humidity of 80%. For these the size frequency distribution is written 


dN 

Id 


M 


_ dNj 
~ dD 
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dD 


( 8 ) 


where, dNi is the number of particles per unit volume per with diameter between D and 
D + dD of the i th species or component. The individual components are taken to be 
log-normally distributed, 12 
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where, Di , and cr; are the modal diameter and the standard deviation, respectively, and 
Ni is the total number density of the t th component. Table 1 provides the parameters Ni, 
Di, <Ti, and m r — irrii for each test aerosol model. The nomenclature M80, C80, T80, and 
U80, refers to the Maritime, Coastal, Tropospheric, and Urban aerosol models 12 at 80% 
relative humidity. Figure 2 provides the spectral variation in r a for the four test aerosol 
models. Note that the spectral variation of the test data is enclosed within the range of 
the candidate aerosol set, for which r a (A)/T 0 (865) « (SfiS/A)^ -2 ). To provide a sense 
of how well the candidate size distribution approximates the test distributions, Figure 3 
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compares the volume size distribution dV/dlnD of three of the four test models (solid 
curves) with that for the Junge power-law distribution (dotted curves). In Figure 3, all 
volume distributions are normalized such that J^° dN = 1. Because of their bimodality, 
M80 and U80 can be better represented by Junge power-law distributions than T80. Table 
2 provides the single-scattering albedo for the test models. Clearly, U80 absorbs much 
more strongly than the others, which we refer to as “weakly” absorbing. 

A test data set was simulated by computing p t (A) — t(A)/9u,(A) for the M80, C80, T80, 
and U80 aerosol models with r a (865) = 0.1, 0.2, and 0.3, and p w (A) with 6° = 0.3 m -1 
(the adopted mean value 3 for Case 1 waters) and C = 0.1, 0.5, and 1.0 mg/m 3 . Note 
that for the test data, the aerosol optical thickness at 412 nm can be as high as 0.80 for 
T80 (Figure 2). The simulated data were created for solar zenith angles (# 0 ) of 0°, 20°, 
40°, and 60°, with the sensor viewing at an angles ( 9 ) of ss 1° and ss 45° with respect 
to nadir and an azimuth ( 4 > ) of 90° relative to the sun. These cover much of the sun- 
viewing geometry available to an ocean color sensor. The combination 9 0 = 0 and 9 ss 1° 
is excluded, as it would be within the sun’s glitter pattern, so there are a total of seven 
sun- viewing geometries considered. 

B. Retrieval results for error-free p t . 

In this subsection, we examine the performance of the algorithm when the TOA 
reflectance p t is free of error. Examples of the retrieval of the spectra of p a + p ra and 
ip w , along with the true values and the recovered values of the parameters, are provided 
in Figure 4 for two different sun-viewing geometries and pseudodata generated using two 
different aerosol models (M80 and U80). The retrievals were effected using the first five 
SeaWiFS bands (412, 443, 490, 510, and 555 nm) in the nonlinear optimization, and 765 
and 865 nm to determine v and r„(865). The results indicate excellent fits between the 
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measured and reconstructed values of p t (A) — p r (A) are obtained for both the strongly 
(U80) and weakly (M80) absorbing aerosols. The retrieved tp w (X) is also excellent, with 
an error at 443 nm less than 0.001 (recall that Gordon and Wang 9 used as a criterion 
for the performance of their SeaWiFS algorithm that the error in tp t at 443 nm should 
be smaller than ±0.002). Retrievals of (7, u>o, and 6° are also excellent. In contrast, the 
retrieval of r a (865) is sometimes very poor, e.g., Figure 4c and 4d. This failure to derive 
acceptable values of r a (865) is completely explained by the fact that the scattering phase 
functions resulting from the Junge power-law distribution do not correspond well to those 
used in the generation of the pseudodata. Figures 5a and 5b compare the phase functions 
[P(0), where 0 is the scattering angle] derived in the retrievals, i.e., computed using the 
retrieved values of is, rrii , and m r , with those used in the generation of the pseudodata 
for Figures 4c and 4d respectively. There are significant differences between the retrieved 
and true phase functions, and for the geometry in this situation (single scattering angle 
0 = 120°) the difference causes the low values of r a (865) obtained by the algorithm. 

Figure 6 provides similar results for pseudodata created using the T80 aerosol model. 
In this case the results are still acceptable, but not as good as those in Figure 4. The degra- 
dation in this case is caused partially by the fact that the Junge power-law distribution 
cannot fit the T80 distribution well (Figure 3), and partially by the fact that the optical 
thickness in the visible can become very large (Figure 2) causing the reflectances to become 
very large (compare Figure 6b with Figures 4c and 4d). In spite of the deficiencies of the 
algorithm, it is remarkable that in Figures 4 and 6, the maximum error in the retrieved 
value of <7 is only slightly over 10%. 

The performance of the algorithm for other geometries, pigment concentrations, and 
aerosol optical thicknesses is summarized in Tables 3-5. In these tables, the quantity being 
retrieved (a> 0 ) C, and b° for Tables 3, 4, and 5, respectively) is averaged over the seven 
sun viewing directions, and its standard deviation (5j) over the seven directions is also 
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presented along with the error in the averaged quantity ( D ). A small value for Sd indicates 
stability in the performance as the sun-viewing geometry is varied. Clearly, it is desirable 
for both Sd and D to be as small as possible. 

The results for u>o at 865 nm (Table 3) show that weakly- and strongly-absorbing 
aerosols are well distinguished; however, the actual derived value of ujq is not as accurate 
as that obtained for the same test data using the Gordon, Du, Zhang 11 algorithm (Table 
2 in Ref. 11). In particular, larger errors are apparent for the weakly-absorbing aerosols; 
however, the retrieval is generally improved as r 0 (865) increases, i.e., as the atmospheric 
contribution to p t increases. 

The retrievals of C (Table 4) are generally excellent, with the largest error being 
17.6%. However, the T80 results are poorer than the others — the largest error in the 
absence of T80 being 10.4%. The error in C increases with increasing r a (865), but not 
dramatically. Other than a few T80 cases, Sd is generally < 5% indicating good stability 
with respect to sun- viewing geometry. Table 5 shows that the retrieval of 6° can usually be 
carried out with an error less than 20% (except for T80) with good sun- viewing geometry 
stability. 

We have not tabulated the performance regarding the retrievals of r a (865), as they are 
significantly poorer than those in Gordon, Du, and Zhang 11 due to the strong differences 
in the Shettle and Fenn 12 and the Junge power-law phase functions (Figure 5). 

The quality of the resulting fits of the optimized and TOA pseudoreflectances is mea- 
sured by the value of Slsq ■ Table 6 provides the mean value of Slsq and its standard 
deviation as over all values of r 0 (865) and C provided in Tables 3-5, as well as over all 
sun- viewing geometries (a total of 63 individual determinations for each model used to 
generate the pseudodata). [The models UU80, U280, and U480 are described later.] The 
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table shows that excellent fits are obtained with relatively low dispersion. As expected, 
the error is somewhat larger for the T80 case. 

The results of this subsection suggest that the water parameters C and b° can be 
retrieved with reasonably good accuracy. The poorer results for u> 0 and, in particular, 
for r a (865) are simply manifestations of the fact that, as argued by Wang and Gordon, 27 
an aerosol model that yields a phase function that correctly approximates the true phase 
function is required to retrieve accurate values, even for weakly-absorbing aerosols. How- 
ever, this is apparently not the case if the goal is only to retrieve the water parameters, 
i.e., a correct representation of the aerosol phase function is not necessary for retrieving 
water parameters. 

C. Retrieval results with calibration error in p t . 

To illustrate the influence of radiometric calibration error on the performance of this 
algorithm, we follow Gordon, Du, and Zhang 11 and add error to pt equal to the uncertainty 
expected after effecting an in-orbit calibration adjustment of the visible bands with respect 
to the band at A/. Gordon 28 has shown that for such a procedure, the residual calibration 
error in all of the bands will have the same sign as that at A/, and the magnitude of the 
error will progressively decrease into the visible, roughly as p r increases. Table 7 provides 
the estimated residual uncertainty (in %) given an uncertainty of 5% in the radiometry 
at A/. We have added (subtracted) the error in Table 7 to (from) the pseudodata used in 
Subsection 3B to produce erroneous reflectances. The algorithm was then operated using 
these reflectances as input pseudodata. Figure 7 provides examples of the retrievals in the 
presence of these errors. It shows that the retrieval of p a -f- p ra is too large (small) if the 
radiometric error is positive (negative). Typically, this is compensated by the opposite 
behavior in tp w . To achieve this the algorithm apparently varies 6° holding C to the 
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nearly-correct value (compare Figures 7a and 7b with Figure 4c, and Figures 7c and 7d 
with Figure 4d). Note that the calibration error will cause an incorrect v to be chosen: 
a positive (negative) error results in a smaller (larger) u. The quality of the fits is only 
marginally degraded from the error-free case. Averaging Slsq for the M80, C80, T80, and 
U80 in Table 6 yields 0.26%. This is increased only to 0.37% for a negative calibration 
error (the worse case of the two). 

Tables 8 and 9 summarize all of the retrievals in the presence of radiometric calibration 
error. They should be compared to Table 4. Although there are some systematic variations 
in C with calibration error, e.g., for C = 0.1 mg/m 3 the retrieved C tends to be larger 
(smaller) for a positive (negative) calibration error, it is not systematic over the range 
of C values or aerosol models. We can conclude from this analysis that the algorithm 
is relatively insensitive to radiometric calibration error of the magnitude expected after 
in-orbit calibration adjustment. It is certainly no more sensitive than the Gordon, Du, 
and Zhang 11 algorithm. 

D. Influence of aerosol vertical structure on retrieval of C 

In the results described in Subsections 3B and 3C, it has been assumed that the 
aerosol is located in a thin layer near the sea surface, as it has been shown that the 
vertical structure of the aerosol concentration is irrelevant as long as the aerosol is weakly 
absorbing. 10 However, when the aerosol is strongly absorbing, this is not the case. 10 ’ 11 In 
the earlier subsections, the vertical structure of the test pseudodata and the candidate 
Junge power-law models were identical. In this subsection, we examine the quality of the 
retrievals when the Urban aerosol is uniformly mixed with air (molecular scattering) from 
the surface to 2 km (U280), from the surface to 4 km (U480), and uniformly mixed with air 
from the surface to the TOA (UU80). In contrast, the radiative transfer computations for 
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the candidate aerosol models still employ a two-layer atmosphere with all aerosol scattering 
in the lower layer and all molecular scattering in the upper layer. 

Figure 8 provides examples of the performance of the algorithm for the U280 and U480 
cases with 0q = 20° and 0 = 1.02°. These should be compared with Figure 4b for which the 
pseudodata and the candidate models had the same vertical distribution of aerosol. Note 
the decreasing quality of the retrievals as the aerosol layer thickens. This is also manifest 
in a decreasing quality of the fit of the computed reflectances to the pseudodata, i.e., 
increasing Slsq> with increased layer thickness (Table 6). Table 10 provides a summary 
comparing the retrieved <7, averaged over sun-viewing geometry, for the U280 and U480 
with that for U80. The UU80 retrievals are not shown because the maximum value of C 
allowed (1.5 mg/m 3 ) was always chosen. Note that as the layer thickens, the retrieved 
values for C increase. This is explained as follows. Increasing the layer thickness causes 
/>a(A) + Pro(A) to become increasing smaller as A decreases (see Figure 4 of Ref. 11). The 
algorithm tries to cope with this decrease in absorption by (1) choosing the largest value of 
rrii available, and (2) increasing C to further reduce p*(A) -p r (A) in the blue relative to the 
red. The fact that the U280 retrievals were realistic is encouraging. It suggests that the 
algorithm as presented could handle aerosol vertical structure with only minor changes. 
The reasoning is as follows. First, as most strongly absorbing aerosols, e.g., wind-blown 
dust, smoke from biomass burning, etc., are likely to be mixed to altitudes of 2-4 km, it is 
clear that the simple vertical structure of the atmosphere used in generating the candidate 
model LUTs (a two-layer model: all aerosol scattering in the lower layer and all molecular 
scattering in the upper layer) must be modified. Second, since the vertical structure in the 
case of nonabsorbing or weakly-absorbing aerosols is irrelevant, using LUTs with aerosol 
mixed higher in the atmosphere will not have a negative impact on the retrieval of C for 
such aerosols (M80, C80, and T80). Third, the U280 retrievals were good even though 
the vertical distribution of the aerosol mixing was in error by 2 km. These observations 
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suggest that if the candidate LUTs were all generated for an aerosol mixed to 2 km, the 
performance in the case of the M80, C80, T80, and U80, examples would be only slightly 
degraded, from that in Table 4, while the performance for U480 would be significantly 
improved over that in Table 10 and that for U280 would be slightly improved. This is the 
strategy that will be used in implementing the algorithm. 

4. Discussion. 

The results in this paper are in agreement with the conclusion of Zhao and Nakajima: 14 
the Junge power-law aerosol size distribution can be usefully applied to the problem of at- 
mospheric correction of ocean color sensors. However, we have extended their algorithm, 
in which only a single value of refractive index (1.50 — iO.Ol) was used for both the Junge 
power-law distribution and the test pseudodata, to allow both m r and rrii to be determined 
in the nonlinear optimization procedure. Simulations show that our algorithm yields ex- 
cellent retrievals of C in atmospheres with either weakly- or strongly-absorbing aerosols, 
but retrieval of the aerosol optical thickness can show large error. The error in r 0 ( 865) 
is traced to significant differences between the phase functions of the optimized Junge 
power-law and the test pseudodata (Figure 5). In contrast, Zhao and Nakajima typically 
obtained excellent r a (550) retrievals (error ~ 10%). This discrepancy is probably due to 
the fact that they used the same refractive index for all of their computations. Hansen 
and Travis 16 show that for a fixed refractive index, the phase functions computed for sig- 
nificantly different size distributions can be very similar as long as the size distributions 
have the same effective radius and effective variance. In our more realistic simulations, in 
which even the individual components of the log-normal size distributions for the test data 
[Eqs. (8) and (9)] have different refractive indices [which also depend on A], this similarity 
is precluded, preventing an accurate estimate of r tt (865). 
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Our results suggest that the detailed shape of the aerosol scattering phase function 
is not relevant to atmospheric correction in Case 1 waters, i.e., we can retrieve excellent 
values of C (Figure 4) with very poor phase functions (Figure 5). Apparently, the most 
important parameters are oj 0 and the spectral variation of optical thickness, r 0 (A)/r 0 (865). 
[Note that the absolute value of r a (865) is directly dependent on the phase function which 
is not relevant.] Gordon 10 has already pointed out that the choice of candidate aerosol 
models in the Gordon and Wang 9 correction algorithm forces a particular u>o-T a (A)/r 0 (865) 
relationship to be applied in the correction. In the present algorithm, no such relationship 
is impressed on the retrievals — v 0 is determined mostly by u and mi, and r a (A)/r a (865) 
is principally determined by u. 

The fact that excellent retrievals of C are obtained even though the phase function of 
the resulting aerosol model is a poor approximation to the true phase function (Figure 5) 
suggests that this algorithm is not limited by the fact that the aerosol optical properties 
[wo, T 0 (A)/r a (865), and P(0)] are computed assuming the aerosol consists of homogeneous 
spheres (Mie theory). Mishchenko and co-workers 29 ’ 30 have shown that in the case of 
randomly oriented spheroids, used as a model for dust-like tropospheric aerosols, the only 
aerosol optical property that is strongly influenced by the particle shape is the phase 
function. Thus, for the purpose of retrieving C, the use of Mie theory should not restrict 
the applicability of the algorithm. In contrast, the retrieved r a (865) depends directly on 
P(Q) and therefore will only be correct if P(Q) is correct. 27 

We have made the implicit assumption that m, and m r are constant for the Junge 
power-law models. This results in only a small variation in u>o with A (Table 2). Unfor- 
tunately, some aerosol species, such as wind-blown dust, are colored and show significant 
variation in m, (and therefore u>o) with wavelength. 31 The algorithm as described above 
cannot handle these aerosols; however, with minor changes it should be able to. Essen- 
tially, all that is needed is to force a prescribed spectral variation in m r and rtii — which, 
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along with u and r 0 (865) determined as in Subsection 3A, fixes all of the aerosol properties 
— and vary only C and b 0 until a minimum in Slsq is reached. If the correct spectral 
variation in m r and m; is chosen, the minimal Slsq in this case should be considerably 
smellier than that obtained using wavelength-independent values of m r and m*. Thus, a 
larger-than-expected minimal Slsq could be employed as a signal that colored aerosols 
may be present, and that a set of variable-m values should be tested to see if Slsq can be 
reduced. 

Our proposed method for dealing with the vertical distribution of strongly-absorbing 
aerosols — placing them in a uniformly mixed layer with a thickness of 2 km in the 
computations used to generate the Junge power-law LUTs — may not be optimum. Un- 
fortunately, other approaches, such as trying to retrieve the layer thickness using LUTs 
generated for several layer thicknesses, 11 may require too much computer memory to be 
useful. There were 72 separate aerosol models already used as the basis for our interpo- 
lation of p a + Pro for arbitrary u, m r , and m*. These require about 500 Mb of random 
access memory to operate the algorithm. In any event, in the final analysis, all atmo- 
spheric correction procedures must compromise in some manner to deal with the problem 
of absorbing-aerosol vertical structure. 

We believe that the principal advantage of the present algorithm over and above 
that proposed by Gordon, Du, and Zhang 11 is the fact that realistic multi-component 
aerosol models are not required for the extraction of ocean properties (although the spectral 
variation of m will be required for colored aerosols). Unfortunately, by abandoning realistic 
models in favor of the simpler Junge power-law, we lose some ability to retrieve meaningful 
aerosol properties, e.g., r a (865). 

We are presently studying ways of making this algorithm more efficient, with the goal 
that it can be applied to the processing of SeaWiFS and MODIS imagery. 
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Figure Captions 


Figure 1. Relationship between £'(765, 865) and v (solid curve) averaged over over all 
combinations of four values of r a (865) (0.1, 0.2, 0.3, and 0.4), six values of m; (0, 0.001, 
0.003, 0.01, 0.03, and 0.04), and two values of m v (1.333, and 1.50), with 0q = 20° and 
0 = 45.9°. Dashed curves represent the upper and lower envelopes of the relationship and 
the dark speckled area represents the standard deviation of the various cases about the 
mean for the given value of v. 

Figure 2. Spectral variation of the aerosol optical thickness for the models used to generate 
the test pseudodata. 

Figure 3. Comparison between the volume size distribution of the aerosol models used to 
generate the test pseudodata (solid curves) and the Junge power-law distribution (dotted 
curves). 

Figure 4. Comparison between the given (curves) and the retrieved (symbols) reflectances 
using the Junge power-law size distribution for 0 = 1.02° (i.e., near nadir): (a) M80 with 
0 O = 20°, (b) U80 with 0 O = 20°, (c) M80 with 9q = 60°, (d) U80 with 9q = 60°. Unprimed 
parameters are the true values, primed parameters are the retrieved values. 

Figure 5. Comparison of the true (solid curve) and retrieved (dashed curve) phase functions 
at 865 nm. The retrieved phase function is used to derived the aerosol optical thickness in 
Figures 4c and 4d: (a) M80, (b) U80. 

Figure 6. Comparison between the given (curves) and the retrieved (symbols) reflectances 
using the Junge power-law size distribution for 0 = 1.02° (i.e., near nadir): (a) T80 with 
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9o — 20°, (b) T80 with 9q — 60°. Unprimed parameters are the true values, primed 
parameters are the retrieved values. 

Figure 7. Comparison between the given (curves) and the retrieved (symbols) reflectances 
using the Junge power-law size distribution for 9 = 1.02° (i.e., near nadir) in the presence 
of the calibration errors in Table 6: (a) M80 with 9 0 = 60° and positive calibration 
error, (b) M80 with 9 0 = 60°, and negative calibration error, (c) U80 with 6 0 = 60° 
and positive calibration error, (d) U80 with 9q = 60°, and negative calibration error. 
Unprimed parameters are the true values, primed parameters are the retrieved values, and 
p c represents the calibration error. 

Figure 8. Comparison between the given (curves) and the retrieved (symbols) reflectances 
using the Junge power-law size distribution for 9 0 = 20° and 9 — 1.02° (i.e., near nadir) 
for cases in which the aerosol is uniformly mixed with air from the surface to an altitude h: 
(a.) U80 aerosol with h = 2 km (U280), (b) U80 aerosol with h = 4 km (U480). Unprimed 
parameters are the true values, primed parameters are the retrieved values. 
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Figure 1. Relationship between £'(765,865) and v (solid curve) 
averaged over over all combinations of four values of r a (865) (0.1, 
0.2, 0.3, and 0.4), six values of m* (0, 0.001, 0.003, 0.01, 0.03, 
and 0.04), and two values of m r (1.333, and 1.50), with d 0 = 20° 
and 6 = 45.9°. Dashed curves represent the upper and lower en- 
velopes of the relationship and the dark speckled area represents 
the standard deviation of the various cases about the mean for 
the given value of v . 




Figure 2. Spectral variation of the aerosol optical thickness for 
the models used to generate the test pseudodata. 
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Figure 3. Comparison between the volume size distribution of 
the aerosol models used to generate the test pseudodata (solid 
curves) and the Junge power-law distribution (dotted curves). 
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Figure 5. Comparison of the true (solid curve) and retrieved (dashed curve) phase functions at 865 
nm. The retrieved phase function is used to derived the aerosol optical thickness in Figures 4c and 
4d: (a) M80, (b) U80. 





(a) (b) 


Figure 6. Comparison between the given (curves) and the retrieved (symbols) reflectances using 
the Junge power-law size distribution for 6 — 1.02° (i.e., near nadir): (a) T80 with 6 0 = 20°, (b) 
T80 with &o = 60°. Unprimed parameters are the true values, primed parameters are the retrieved 
values. 
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Figure 7. Comparison between the given (curves) and the retrieved (symbols) reflectances using the 
Junge power-law size distribution for 0 = 1.02° (i.e., near nadir) in the presence of the calibration 
errors in Table 6: (a) M80 with 60 = 60° and positive calibration error, (b) M80 with 6 q = 60°, 
and negative calibration error, (c) U80 with 9 0 = 60° and positive calibration error, (d) U80 
with 9q = 60°, and negative calibration error. Unprimed parameters are the true values, primed 
parameters are the retrieved values, and p c represents the calibration error. 
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00 


Figure 8. Comparison between the given (curves) and the retrieved (symbols) reflectances using 
the Junge power-law size distribution for d 0 = 20° and 6 = 1.02° (i.e., near nadir) for cases in 
thich the aerosol is uniformly mixed with air from the surface to an altitude h: (a) U80 aerosol 
with h = 2 km (U280), (b) U80 aerosol with h = 4 km (U480). Unprimed parameters are the true 
values, primed parameters are the retrieved values. 





Table 1: Characteristics of the test aerosol models used in the study. 


Aerosol Model 

Size Distribution 

Refractive Index 


N, 

D, 

<Ti 

412 nm 

865 nm 

M80 

0.990000 

0.06548 

0.35 

1.446 — t'3.309E— 3 

1.436-16. 107E-3 


0.010000 

0.63600 

0.40 

1.359 — t'5.165E— 9 

1.348 — tl.38lE— 6 

C80 

0.995000 

0.06548 

0.35 

1.446— i'3.309E— 3 

1.436-i6. 107E-3 


0.005000 

0.63600 

0.40 

1.359— t5.165E— 9 

1.348— tl.38lE— 6 

T80 

1.000000 

0.06548 

0.35 

1.446 — t'3.309E— 3 

1.436 — {6.107E- 3 

U80 

0.999875 

0.07028 

0.35 

1.423— i'3.473E— 2 

1.414— i3.412E— 2 


0.000125 

1.16200 

0.40 

1.415 — t'3.15lE— 2 

1.406— t'3.095E— 2 


Table 2: Values of single scattering albedo at A = 412 and 865 nm 
for the test aerosol models used in the study. 


Aerosol Model 

A (nm) 

Vo 

M80 

412 

0.9923870 


865 

0.9934230 

C80 

412 

0.9883920 


865 

0.9884390 

T80 

412 

0.9758390 


865 

0.9528370 

U80 

412 

0.7823030 


865 

0.7480590 





Table 3 


RETRIEVAL OF AERO SOU s SINGLE SCATTERING ALBEDO 

Tabulated are mean values of retrieved G) 0 for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80, C80 , T80, U80). 

Also, provided are the standard deviations over viewing geometries divided by the mean ( S d ) 
as well as the deviations from given parameters ( D ). 



C: mg/m^ 

0.100 

0.500 

1.000 



result 

s d 

D 

result 

s d 

D 

result 


D 

M80 

x a (865) = 0.10 

0.938 

2.5% 

5.50% 

0.968 

1.9% 

2.40% 

0.972 

2.2% 

2.12% 

C0 0 = 0.993 

t a (865) = 0.20 

0.954 

2.2% 

3.93% 

0.971 

2.4% 

2.22% 

0.974 

2.5% 

1.91% 


t a (865) = 0.30 

0.963 

1.8% 

3.02% 

0.976 

1.9% 

i .22% 

0.979 

2.0% 

1.41% 



C: mg/m^ 

0.100 

0.500 

1.000 


result 

s d 

D 

result 

s d 

D 

result 

s d 

D 

C80 

t a (865) = 0.10 

0.943 

2.3% 

4.56% 

0.970 

1.7% 

3.20% 

0.972 

1.7% 

1.62% 

O> 0 =0.988 

t a (865) = 0.20 

0.960 

2.1 % 

2.83% 

0.977 

1.7% 

i.i i% 

0.979 

1.8% 

0.91% 


r a [865) = 0.30 

0.969 

1.9% 

1 .92% 

0.981 

1.6% 

0.71% 

0.983 

1.5% 

0.51% 



C; mg/m^ 

0.100 


1.000 



result 

s d 

D 


m 


result 

S d 

D 

T80 

x a (865) = 0.10 









8.50% 

to,, = 0.953 

X a (865) = 0.20 



IH 



6.51% 





x a (865) = 0.30 

0.902 

1.6% 



1.0% 

4.09% 



3.67% 



C; mg/m^ 

O.lOO 

0.500 

1.000 





D 

result 

s d 

D 

result 


D 


x a (865) =0.10 



0.54% 

0.745 

4.3% 

1.40% 

0.745 




x a (865} = 0.20 



2.94% 

0.769 

0.8% 

2.81% 

0.770 

1.3% 

2.94% 


x a {865) = 0.30 

0.788 


5.35% 

0.777 

1.2% 

3.88% 

0.773 

1.5% 

3.34% 
















Table 4 


RETRIEVAL OF OCEAN'S PIGMENT CONCENTRATION 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80, C80, T80, U80). 

Also , provided are the standard deviations over viewing geometries divided by the mean ( S d ) 
as well as the deviations from given parameters ( D ). 



C: mg/ m 3 

0.100 

0.500 

1.000 


result 

s d 

D 

result 

s d 

D 

result 

s d 

D 

M80 

t a (865) = 0.10 

0.096 

1.8% 

4.00% 

0.488 

1.8% 

2.40% 

1.012 

1.1% 

1 .20% 

® 0 = 0.993 

x a (865) = 0.20 

0.094 

2.9% 

6.00% 

0.479 

3.7% 

4.20% 

1.014 

2.3% 

1 .40% 


t a (865) = 0.30 

0.093 

3.2% 

7.00% 

0.473 

4.6% 

5.40% 

1.010 

3.5% 

1.00% 



C; mg/m 3 

0.100 

0.500 

1.000 



result 

S d 

D 

result 

s d 

D 

result 

S d 

D 

C80 

x a (865) = 0.10 

0.096 

2.1% 

4.00% 

0.484 

2.1% 

3.20% 

1.001 

1.9% 

0.10% 

(D u =0.988 


0.094 

3.2% 

6.00% 

0.469 

1.7% 

6.20% 

0.984 

4.3% 

1 .60% 


t a (865) = 0.30 

0.093 

2.8% 

7.00% 

0.459 

4.9% 

8.20% 

0.959 

7.3% 

4.10% 



C; mg/m 3 

0.100 

0.500 

1.000 





D 

result 

S d 

D 

result 

S d 

D 

T80 

t a (865) = 0.10 

HI 



0.466 

4.0% 

6.80% 

1.021 

6.1% 

2.10% 

CO u = 0.953 

t a (865) = 0.20 

1 



0.425 

9.9% 

15.0% 

0.941 

7.0% 

5.90% 


t a (865) = 0.30 

1 



0.412 

6.8% 

17.6% 

0.874 

6.7% 

12.6% 



C: mg/m 3 

0.100 

0.500 

1.000 



result 

S d 

D 

result 

s d 

D 

result 

s d 

D 

U80 

t a (865) = 0.10 

0.100 




0.7% 


1.028 

1.9% 


= 0.748 

x a (865) = 0.20 




0.526 




3.2% 



x a (865) = 0.30 



10.0% 

0.552 






























Table 5 


RETRIEVAL OF WATER SCATTERING PARAMETER 

Tabulated are mean values of retrieved b ° for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80, C80, T80, U80). 

Also, provided are the standard deviations over viewing geometries divided by the mean ( ) 

as well as the deviations from given parameters ( D ). 



C: mg/m 3 

0.100 

0.500 

1.000 



result 


D 

result 

s d 

D 

result 

c 

D 

M80 

t a (865) = 0.10 

0.294 

1.1% 

2.00% 

0.274 

4.1% 

8.67 % 

0.273 

4.4% 

3.20% 

(0 0 = 0.993 

t a (865) = 0.20 

0.289 

1.7% 

3.67% 

0.258 

7.3% 

14.% 

0.255 

8.2% 

15.0% 


x a (865) =0.30 

0.286 

1.8% 

4.67% 

0.249 

8.2% 

17.0% 

0.246 

9.4% 

18.0% 



Cl mg/m 3 

0.100 

0.500 

1.000 



result 

S d 

D 

result 

s d 

D 

result 

El 

D 

C80 

x a (865) = 0.10 



2.33% 

0.273 

4.9% 

3.20% 



9.67% 

to,, =0.988 

x a (865) = 0.20 



3.67% 

0.255 

8.0% 

15.0% 



16.7% 


x a (865) = 0.30 

1 


4.33% 

0.247 

7.4% 

17.7% 

0.239 

8.7% 

20.3% 



C: mg/m 3 

0.100 

0.500 

1.000 



result 

S d 

D 


S d 

D 


wa 

D 

T80 

t a (865) = 0.10 





18.% 

21.0% 




a> 0 = 0.953 

x a (865) =0.20 





23.% 

33.3% 


23.79 

3 36.% 


x a C865) = 0.30 





20.% 

33.0% 


25.% 




C; mg/m 3 

0.100 

0.500 

1.000 



result 

wm 

D 



D 


E9 

D 

U80 

x a (865) = 0.10 

0.297 





0.67% 


1.7% 

0.67% 

to,, = 0.748 

x a (865) = 0.20 

0.302 




1.9% 

5.33% 


2.4% 

6.33% 


x a (865) = 0.30 

0.311 







8.1% 

18.0% 

































Table 6 : Mean Slsq and its standard deviation as (both in %), 
over all values of r 0 (865), C , and sun-viewing geometries. 


Model 

Slsq 

<rs 

M80 

0.27 

22 

C80 

0.28 

19 

T80 

0.36 

23 

U80 

0.15 

16 

U280 

1.57 

27 

U480 

5.58 

27 

UU80 

133.48 

68 




Table 7: Values of the residual radiometric calibration uncertainty 
after effecting an in-orbit calibration adjustment. 28 


Ai 

(nm) 

% Uncertainty 

412 

0.3 

443 

0.5 

490 

0.8 

520 

1.0 

550 

1.5 

670 

2.0 

765 

3.0 

865 

5.0 



Table 8 


RETRIEVAL OF OCEAN'S PIGMENT CONCENTRATION 

with negative calibration error 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80, C80, T80, U80). 

Also , provided are the standard deviations over viewing geometries divided by the mean ( S d ) 
as well as the deviations from given parameters ( D ). 



C; mg/m^ 

0.100 

0.500 

1.000 


result 


D 

result 

s d 

D 

result 

s d 

D 

M80 

t a (865) =0.10 

0.094 

2.3% 

6.00% 

0.479 

2.3% 

4.20% 

1.004 

0.8% 

0.40% 

(0 0 = 0.993 

x a (865) = 0.20 

0.091 

4.0% 

9.00% 

0.463 

4.6% 

7.40% 

1.005 

2.2% 

0.50% 


t a (865) = 0.30 

0.090 

4.7% 

3.00% 

0.452 

6.1% 

9.60% 

0.998 

4.0% 

0.20% 



C; mg/m^ 

0.100 

0.500 

1.000 



result 

S d 

D 


s d 

D 

result 


D 

C80 



2.6% 



2.5% 

4.80% 

0.997 



10 0 = 0.988 

x a (865) = 0.20 

ihi 

4.7% 



5.2% 

9.00% 


2.2% 

1.70% 


x a (865) = 0.30 

m 

5.3% 


0.438 

6.5% 

12.4% 

0.964 

4.4% 

3.60% 



C; mg/m^ 

0.100 

0.500 




result 

s d 

D 

result 

S d 

D 


s d 

D 

T80 

x a (865) = 0.10 

0.086 

10.% 

5.00% 

0.437 

9.5% 

12.6 % 


14.% 

3.80% 

cd 0 = 0.953 

x a (865) = 0.20 

0.078 

14.% 

22.0% 

0.377 

19.% 

24.6% 

0.836 

3.4% 

16.4% 


X a (865) = 0.30 

0.081 

12.% 

19.0% 

0.337 

15.% 

32.6% 

0.640 

26.% 

36.0% 



C: mg/m^ 

0.100 

0.500 

1.000 





D 



D 

result 

s d 

D 

1 

x a (865) = 0.10 


1.5% 


w 


0.80% 




1 

x a (865) = 0.20 

0.099 


ESI 


1.9% 


1.033 


B 


x a (865) =0.30 

0.104 




4.7% 


1.056 

4.1% 

5.60% 



















Table 9 


RETRIEVAL OF OCEAN'S PIGMENT CONCENTRATION 
with positive calibration error 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80, C80, T80, U80). 

Also , provided are the standard deviations over viewing geometries divided by the mean ( ) 

as well as the deviations from given parameters ( D ). 



C; mg/m 3 

O.JOO 

0.500 

1.000 


result 

s d 

D 

result 


D 

result 

S d 

D 

M80 

t a (865) = 0.10 

0.098 

1.6% 

2.00% 

0.494 

2.7% 

1.20% 

1.011 

4.2% 

1.10% 

C0 0 = 0.993 

T a (865) = 0.20 

0.097 

2.2% 

3.00% 

0.487 

5.0% 

2.60% 

0.999 

8.5% 

0.10% 


x a (865) =0.30 

0.097 

2.2% 

3.00% 

0.485 

6.2% 

3.00% 

0.991 

11.% 

0.90% 



C: mg/m 3 

0.100 

0.500 

1.000 



result 

s d 

D 

result 

s d 

D 

result 

s d 

D 

C80 

t a (865) = 0.10 

0.098 

1.6% 

2.00% 

0.490 

3.6% 

2.00% 

0.995 

6.3% 

0.50% 

co 0 =0.988 

x a (865) = 0.20 

0.097 

2.2% 

3.00% 

0.478 

7.1% 

4.40% 

0.966 

12.% 

3.40% 


x a (865) = 0.30 

0.097 

2.6% 

3.00% 

0.471 

10.% 

5.80% 

0.940 

17.% 

6.00% 



C: mg/m 3 

0.100 

0.500 

1.000 



result 

s d 

D 

result 


D 

result 



T80 


0.095 

4.7% 

5.00% 


2.4% 

2.60% 




(0 0 = 0.953 

x a (865) = 0.20 

0.093 

1.3% 

7.00% 


4.0% 

6.00% 

0.935 




x a (865) = 0.30 

0.096 

4.4% 

4.0% 


2.1% 

6.20% 

0.943 

2.7% 

5.70% 



C; mg/m 3 

O.JOO 



1.000 




resu/t 

s d 

D 


S d 

D 

result 

S d 

D 

U80 

x a (865) = 0.10 



IHI 

V 

1.1% 

3.80% 




C0 () = 0.748 

x a (865) = 0.20 


2.8% 

B 


3.8% 

8.60% 





x a (865) = 0.30 

0.116 



1 

6.4% 

14.6% 


3.6% 



















Table 10 


RETRIEVAL OF OCEAN'S PIGMENT CONCENTRATION 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80, C80, T80 , U80). 

Also, provided are the standard deviations over viewing geometries divided by the mean ( S d ) 
as well as the deviations from given parameters ( D ). 



C: mg/m^ 

0.100 

0.500 

1.000 



result 

s d 

D 

result 

s d 

D 

result 

s d 

D 

U80 

x a (865) = 0.10 

0.100 

0.8% 

0.00% 

0.507 

0.7% 

1.40% 

1.028 

1.9% 

2.80% 

G) 0 = 0.748 

t a (865) = 0.20 

0.104 

1.4% 

4.00% 

0.526 

2.7% 

5.20% 

1.053 

3.2% 

5.30% 


t (865) =0.30 
a 

0.110 

4.3% 

10.0% 

0.552 

5.8% 

10.4% 

1.070 

4.0% 

7.00% 










